Self-consistent random phase approximation in the Schiitte-Da-Providencia 

fermion-boson model 
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The Self-consistent random phase approximation (SCRPA) is applied to the exactly solvable 
model with fermion boson coupling proposed by Schiitte and Da-Providencia. Very encouraging 
results in comparison with the exact solution of the model for various observables are obtained. The 
transition from the normal phase to the phase with a spontaneously broken symmetry is carefully 
investigated. The strong reduction of the variance in SCRPA vs HF is pointed out. 
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During the last decade the so-called Self-Consistent 
version of the Random Phase Approximation (SCRPA) 
has seen very encouraging successes in a number of non 
trivial model cases (see for example Q for a detailed de- 
scription of the method and Q for the application of 
SCRPA to the many level pairing model). In spite of 
these performances of the theory, there are remaining 
problems. In first place this concerns situations with 
spontaneously broken symmetries. Such situations were 
treated in 0, 0, Q. Whereas in the Lipkin model 3] 
the symmetry broken ('deformed') phase caused no prob- 
lem because the broken symmetry is discrete (parity), 
in the other two cases 0, Q , with a continuous broken 
symmetry, problems appeared with the low-lying mode 
known to be exactly at zero energy in the standard HF- 
RPA approach (the spurious or Goldstonc mode) . In 
the two cases cited 0, Q the low-lying mode does not 
appear at zero energy in SCRPA because the RPA oper- 
ator does not contain the symmetry operator as a limit 
case. Indeed, e.g. the number operator in quasiparticle 
(BCS) representation contains a purely hermitian piece 
atcKfe which can not be incorporated in the RPA operator 
which by definition is non hermitian. The same situation 
is present in the Schutte-Da-Providenica boson-fermion 
model where the symmetry operator contains the bo- 
son and fermion number operators. The violation of the 
Goldstone theorem signifies that the Ward identities and 
conservations laws are not respected. Though this viola- 
tion seems relatively mild and to go away in macroscopic 
systems (the hermitian pieces becoming of zero weight), 
the situation remains annoying for finite systems. 

In this note which can be considered as a sequel of |4| 
we want again to investigate the Schutte-Da-Providenica 
model 

H = n + atfb + G(T+tf +r~b) (1) 
with b\ b ideal boson operators and, n — X)«=i a Qi a 0ii 



t° = \{n — n), where the a', a are fermion operators. 
In analogy to the work in Q we will introduce the more 
general RPA operator 

Ql = X v t + -YJ- +KB^ -iXvB+U v fi ] '-V v p/3, v = 1,2, 3 
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The operators T ± , T° are obtained from r , t° by writ- 
ing the latter ones in the "deformed" basis 
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The bosons operators and B are obtained from the 
original ones by a shift transformation B — > b — a \ where 
a is a c-number characterizing the appearance of the 
Bose condensate. The introduction of the boson pair 
operators f3*/3' is motivated by the fact that otherwise 
there exists a certain dissymmetry between fermions and 
bosons, the fermions being in any case bilinear whereas 
the bosons are otherwise only contained to linear order 
in ©. Also the symmetry operator P = tfb — h contains 
the bosons quadratically and the extended ansatz @ 
may therefore show improved behavior with respect to 
the Goldstone mode. The formalism goes exactly in the 
same way as in 0, 0, 01 using the equation of motion 
method 



([SQ,[H',Ql]]) = n v {[6Q,Ql\) 



(5) 
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to determine the amplitudes in © . As in Q , in order to 
fix the value L = (P), we use in (J3J the cranked Hamil- 
tonian H' = H — fiP in the symmetry broken phase, oth- 
erwise H' = H . The mean-field amplitudes u, v and a 
are readily obtained from a minimization of the ground- 
state energy, leading to {[H',t + ]) = {[H',B^]) = 0. The 
amplitudes in form a complete orthonormal set when 
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calculated from |[5J. Then J5J can be inverted and with 
the usual condition for the RPA ground-state 

Q„|RPA) =0, v= 1, 2, 3. (6) 

all expectation values appearing in JQ| like for example 
(t + B) and (B^B) can directly be expressed in 
terms of the RPA amplitudes. The only unknown quan- 
tity at this point preventing a fully self-consistent solu- 
tion of the SCRPA equations Eq. (J5J is the expectation 
value (T }. However, in analogyto our previous study 
for the two-level pairing model [lj this quantity can be 
expressed as an expression in the operators T + and T~ 
up to any order in a fast converging series according to 

N 1 1 TO 

T° = + —T+T- + — -T+ T +.... (7) 

2 N N 2 (N-l) y ' 

With this relation the SCRPA equations are completely 
closed and we can proceed to the numerical solution. We 
notice that with respect to [4J the SCRPA Eq. 101 is a 
6x6 dimensional problem whereas before it was 4x4. 

We now come to the presentation and discussion of 
the results. In what follows, SCRPA(6) refers to the 
SCRPA method with RPA excitation operator quadratic 
in the bosons operators i.e. Eq. J2J), while SCRPA(4) 
refers to the same method but with RPA excitation op- 
erator linear in the bosons operators i.e. U v = V v = 0. 
In the following we also use the set of parameters a = 3, 
N = 30 as in jj] , for which in this model the phase tran- 
sition point is localized at x — 1.0, where x — G\J N/a. 
SCRPA always shows a clear superiority over standard 
RPA, though, besides for some quantities, the differences 
are not very pronounced. Concerning the ground-state 
energy we do not give results but we only notice that we 
arrive practically at the same interpretations as in 
However, in order to test the accuracy of our approach 
it is instructive to calculate the differences of energies of 
the ground-state band with L values just one unit away 
from the absolute ground-state. One such quantity is the 
"chemical potential" which should be identified with the 
Lagrange parameter used for restoring the symmetry 

H=l{E° L+1 -E L _ 1 ). (8) 

In Table H] we show n when we calculate separately El ±1 
(in standard RPA and SCRPA) and then take the differ- 
ence. We also give in Table Q] the L values which corre- 
spond for a given x value to the absolute ground-state. 
In Table[I]we see a strong improvement of SCRPA(6) and 
SCRPA(4) over standard RPA and the high quality of the 
results in comparison with the exact values in the region 
around the phase transition point. Also SCRPA(6) is 
still improved over SCRPA (4). We could also have taken 
the n values found from adjusting the correct L = (P) 
values in the standard RPA and SCRPA calculations; 
we have checked numerically that the results are always 
practically identical. 



X 


L 


fi( exact) 


A*(6) 


M(4) 


At(RPA) 


1.1 


-1 


-0.0131 


-0.0485 


-0.0485 


0.1160 


1.4 


-3 


0.0325 


0.0348 


0.0405 


0.1303 


1.8 


-3 


-0.0341 


-0.0315 


-0.0329 


0.0060 


2.2 





0.0320 


0.0300 


0.0240 


0.0345 


2.6 


4 


0.0305 


0.0285 


0.0221 


0.0207 


3.0 


9 


0.0204 


0.0195 


0.0135 


0.0071 



TABLE I: Chemical potential : exact, SCRPA(6), SCRPA(4) 
and standard RPA; L values, in the deformed region for dif- 
ferent values of the interaction strength x. 

Two other interesting quantities to be calculated 
within the SCRPA formalism and closely related to the 
chemical potential are the energy differences of the ab- 
solute ground-state with its " left" and " right" neighbors 
just one unit in L away 

AE ±1 = ± (E° L±1 El) . (9) 

These quantities are interesting because, as we will ex- 
plain below, they should be closely related to the lowest 
RPA eigenvalue in the symmetry broken phase. Be- 
cause we obtain similar interpretation and conclusions 
for both AE±i, we will present and discuss only the re- 
sult for AE_i. In Fig. ^ we see a very good agreement 
with the exact results of both SCRPA(6) and SCRPA(4). 
However, we note that in this quantity no clear superi- 
ority of SCRPA(6) over SCRPA(4) can be detected, the 
results being at times in favor of the one or the other. 
The good quality of the results for AE±i shows that the 
SCRPA method is able to reproduce the full spectrum. 
We also should notice that the smallness of AE±i means 
that two neighboring ground-states with L and L ± 1, 
respectively, are almost degenerate which indicates that 
the system is in the phase of spontaneously broken sym- 
metry. Furthermore, one can check that in the large N 
limit AE±i tends to zero. The zero eigenvalue (Gold- 
stone mode) which is one of the solutions of the standard 
RPA in the deformed region corresponds to this degener- 
acy of the neighboring ground-state energies in the large 
N limit. 

Let us now discuss the eigenvalues of RPA and SCRPA 
matrices. As it is well known 0,0, in standard RPA the 
lowest eigenvalue corresponds in the symmetry broken 
(deformed) region to the spurious mode Qi — 0, whereas 
the second eigenvalue gives the excitation of the intrinsic 
system. Before coming to this point we should mention 
again that the RPA eigenvalues in the deformed region 
are calculated with the "intrinsic" Hamiltonian H' = H— 
IjlP. Therefore, when the symmetry is restored due to the 
appearance of the Goldstone mode, the RPA eigenvalues 
give the excitation energies of the system. The results for 
the mode VL2 are not shown in a figure because we obtain 
the same interpretation as already given in Q, however, 
with still improved results from SCRPA(6). 

Concerning the low-lying eigenvalue of SCRPA which 
in standard RPA corresponds to the zero-energy eigen- 
value (Goldstone or spurious mode) in the "deformed" 
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FIG. 1: Comparison between the exact, SCRPA(6) and 
SCRPA(4) results for the excitation energy A£_i = E° L — 
E ( [_ 1 in the deformed region. 



FIG. 2: The standard RPA, SCRPA(4) and SCRPA(6) results 
for the spurious mode fii compared with the exact energy of 
the excitation AE—i. 



region, we present in Fig. [21 a comparison between stan- 
dard RPA, SCRPA(4) and SCRPA(6) with exact results. 
In the "spherical" phase we notice that the eigenvalue Oi 
is identified with the exact "intraband" excitation Ai5_i. 
Furthermore, we see the important improvement of the 
SCRPA results in both cases SCRPA(4) and SCRPA(6) 
with respect to standard RPA result. After the phase 
transition A£Li remains finite but very small, slowly de- 
creasing for increasing x, while the lowest eigenvalue in 
standard RPA corresponds to the spurious mode fii = 0. 
Concerning the low-lying eigenvalue in SCRPA calcula- 
tion we see that SCRPA(6) improves the result with re- 
spect to SCRPA(4) but it is still quite far from the ex- 
act result. We can notice that the Goldstone theorem is 
not correctly fulfilled in this case. Therefore, the prob- 
lem of the identification of the low-lying eigenvalue is not 
yet solved in SCRPA in spite of the introduction of the 
quadratic boson terms in Q. 

From the presence of the quadratic boson terms in the 
RPA excitation operator the SCRPA method produces 
an supplementary eigenvalue which is noted O3. Let us 
now discuss the results for this mode, i.e. the third eigen- 
value RPA of the SCRPA(6) which is presented in Fig.|31 
In the spherical region we notice that this mode is iden- 
tified with very good accuracy to the exact "intraband" 
excitation A_E + 2- However, in the deformed region (not 
shown), the result obtained for this mode cannot be iden- 
tified with one of the exact excitations of the system. At 
present we do not have an explanation of this fact. It 
is likely related to the failure of the Goldstone theorem 
mentioned above. See also further discussion of this point 
at the end of this note. The difficulty may be of the same 
origin as with the tti mode. 

A quantity which is particularly sensitive to the correct 
treatment of correlations in the ground-state is the mean 
boson number (not shown in Ref. Q). This expectation 
value can be obtained in terms of the RPA amplitudes 



— AE +2 (exact) 



FIG. 3: The SCRPA(6) results of the energy of the excitation 
mode Q3 compared with the exact energy of the excitation 
AE+2 in the spherical region. 



according to 



N b = (tfb) = (B^B) + a 2 



(10) 



where (B^B) is given, in SCRPA(6), by (B^ B) = /Ltf + 
f4 + f4. In Fig.gjwe show the results of the SCRPA(6), 
SCRPA (4), standard RPA and mean field methods for 
this quantity. Again with SCRPA (6) one notices a sig- 
nificant improvement over standard RPA and HF method 
for which the agreement with the exact result is not sat- 
isfying. Also, we note in Fig. 0] that the SCRPA(6) im- 
proves slightly the result over the SCRPA (4) specially in 
the deformed region. Furthermore, we note that in stan- 
dard RPA method, because we have a Goldstone mode in 
the deformed region, we cannot calculate this quantity. 

Another quantity which is very interesting to investi- 
gate in the SCRPA method is the variance of the sym- 
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metry operator P given by 

AP 2 = (P 2 ) - (P) 2 . (11) 

In Fig. we present the results corresponding to this 
quantity calculated with SCRPA and HF methods. This 
is a new result which was not elaborated in Q. We see 
that the variance is strongly reduced compared to HF 
values. We, however, see that AP even in the SCRPA 
acquires sizable non vanishing values. This simply means 
that the symmetry P, broken at the level of the mean 
field theory, is not completely restored. Furthermore, we 
do not present the standard RPA results concerning this 
quantity in Fig. [S] because the RPA amplitudes originat- 
ing from the Goldstone mode are divergent. This con- 
stitutes the same kind of situation as that for the boson 
number in the ground-state calculation for which we also 
have not given the standard RPA results in deformed 
region. 



physical spectrum. However, in spite of this somewhat 
disappointing result, the introduction of the extra terms 
allowed to reproduce very well a further excited state 



5 



2 



1 1 1 1 1 ' 






— SCRPA 




HF 


-_-/ 1,1,1,1 


, 1,1,1,1,1, 



1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 



X 



0.4 



0.2 



-0.4 





1 1 1 1 1 1 , 

SCRPA(6) 


1 1 1 




■—■ SCRPA(4) 






RPA 






HF 









lili 



I I I I I I I I I I I I I 

0.5 1 1.5 2 2.5 3 

X 



FIG. 4: The difference [N b - N£™ ct ] calculated with 
SCRPA(6), SCRPA(4), RPA and mean field methods as func- 
tion of the interaction strength x. 

In conclusion we reconsidered the work of Bertrand et 
al. 3 who treated the schematic Schutte-Da-Providencia 
model for interacting bosons and fermions within the 
SCRPA scheme. In |4j the RPA operator consisted only 
out of one boson and fermions pairs. Here we extended 
this configuration space and included in addition bosons 
pairs. One of the motivations to do this was to see 
whether the Goldstone theorem which was quite strongly 
violated in 0| is improved. It was found that the low- 
lying mode in the "deformed" zone, i.e. symmetry bro- 
ken region, indeed is lowered by ~ 30% when boson pair 
terms are added to the RPA operator. However, with 
respect to the first physical state, the position of the 
spurious mode is still too high and one therefore can not 
say that it decouples to a good approximation from the 



FIG. 5: Variance as function of the interaction strength x. 



of the spectrum in the symmetry conserved phase and 
additionally the quantities which had already been cal- 
culated in U without the boson pair operators are still 
improved. We also calculated further quantities as the 
number of bosons in the ground-state and the fluctua- 
tion of the symmetry operator. For instance the latter 
becomes strongly reduced with respect to its mean field 
value. However, also for this quantity a substantial non 
vanishing value remains. All in all we can say that the 
inclusion of the extra two boson piece to the RPA oper- 
ator allowed to calculate one more state in the spherical 
region, improve existing results from calculations with- 
out these terms, and lower the spurious state. However, 
with respect to the latter feature no real breakthrough 
could be observed and further ideas are needed to sub- 
stantially improve the situation in the symmetry broken 
region whereas the " spherical" region seems to be well un- 
der control. In this respect this may be similar with other 
approaches treating correlations beyond mean field like 
coupled cluster theory, Jastrow, correlated basis function 
etc. 
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